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Abstract 

We study first the Hamiltonian operator H corresponding to the Fock-Weyl exten- 
sion of the Dirac equation to gravitation. When searching for stationary solutions 
to this equation, in a static metric, we show that just one invariant Hermitian 
product appears natural. In the case of a space-isotropic metric, H is Hermitian 
for that product. Then we investigate the asymptotic post-Newtonian approxima- 
tion of the stationary Schrodinger equation associated with H, for a slow particle 
in a weak-field static metric. We rewrite the expanded equations as one equation 
for a two-component spinor field. This equation contains just the non-relativistic 
Schrodinger equation in the gravity potential, plus correction terms. Those "cor- 
rection" terms are of the same order in the small parameter as the "main" terms, 
but are numerically negligible in the case of ultra-cold neutrons in the Earth's 
gravity. 



Key words: Dirac equation, gravitation, Hamiltonian operator, non-relativistic limit, 
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1 Introduction 

Quantum mechanics in a gravitational field is increasingly becoming an experimentally 
relevant subject Due to the weakness of the gravity interaction, the effects of gravity 
are more clearly seen in experiments with neutral particles: neutrons, and also atoms 
or molecules. The quantum-mechanical phase shift of neutrons due to their interaction 
with the Earth's gravitational field has been measured since a long time, thanks to 
interferometric experiments [2]. {Similar effects, due to a non-inertial (rotational) mo- 
tion or to the Earth's gravity acceleration, have been observed in atom interferometry 
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oil].} More recently, the quantization of the energy levels of ultra-cold neutrons in 
this same terrestrial gravity field has been revealed by measuring their transmission 
through a horizontal slit • New, more precise experiments are being foreseen in or- 
der to accurately determine the energy levels [7] . Up to now, the theoretical analyses of 
these experiments have been done in the framework of the non-relativistic Schrodinger 
equation in the Newtonian gravity potential IHl El HHI • This is justified by the small- 
ness of the velocity of the neutrons used in these experiments (5 m/s to 10 m/s in the 
transmission measurements 010], ca. 3 x 10^ m/s in the interfero metric experiment 
j2j), and by the weakness of the gravitational field of the Earth. 

However, gravity is currently described by relativistic theories like general relativ- 
ity (GR). In the framework of GR and other theories based on Einstein's equivalence 
principle, the wave equations of relativistic quantum mechanics, i.e. essentially the 
Klein-Gordon equation and the Dirac equation, are adapted to gravitation by rewriting 
them in a generally-covariant way. For the Klein-Gordon equation, this "covariantiza- 
tion" is ambiguous ^n], due to the fact that covariant derivatives do not commute. In 
the case of the Dirac equation with the spinor transformation, this procedure leads to 
the equation independently proposed by Fock |12j and by Weyl [12], hereafter named 
the Dirac-Fock-Weyl (DFW) equation. On the other hand, an alternative gravitational 
Dirac equation has been recently proposed (for the case of a static metric) ^1]. It is 
based on a direct derivation from the principles of wave mechanics, instead of using the 
equivalence principle. The new equation definitely violates the latter principle, for it 
does not reduce to the flat-space-time Dirac equation in a local "freely falling" frame 

Since neutrons are particles with spin 1/2, one may hope that, at least in the 
absence of an external electromagnetic field, their behaviour in the gravitational field 
should be correctly described by one of these two gravitational Dirac equations — and, 
of course, it would be interesting to know which one of the two. The aim of this paper 
is to derive, for the standard, DFW equation, a Schrodinger-type equation which will 
allow to compute the stationary energy levels. (The corresponding equation for the 
alternative equation |14J will be derived in a forthcoming work.) The equation to be 
derived exhibits the first correction terms with respect to the non-relativistic stationary 
Schrodinger equation in the Newtonian gravity potential. This does not appear to have 
been done before, although some amount of work has been devoted to studying the 
weak-field and/or non-relativistic limit of the DFW equation (see, among others, Refs. 
|15| ITB] IT7| I18j). These works did result in proposals for an approximate Hamiltonian 
operator, hence one should be able to find easily the approximate equation for the sta- 
tionary energy levels by using those works. However, the approximation scheme was not 
made explicit in these works, so that, in the approximate equations, one cannot easily 
assign orders to the different terms, with respect to a small parameter — as is neces- 
sary to ensure that the first corrections with respect to the non-relativistic Schrodinger 
equation are consistently evaluated. In addition, these works were concerned with the 
Hamiltonian for the four-dimensional complex wave function, whereas we will verify in 
subsect. I2.5l that (not unexpectedly) the solutions of the stationary DFW equation have 
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only two independent complex components. Finally, each of the Hamiltonians which 
were used in Refs. [T31ITEllT7| ll8j. had been obtained from the starting Hamiltonian of 
the DFW equation by using a non- unitary transformation [T^ ITTj. In the absence 
of another interpretation, this procedure might be seen as coordinate-dependent [TU] . 

The paper begins (Sect. ^ with a discussion of the Hamiltonian operator for the 
DFW equation. After noting the invariance of this operator, we investigate the question 
of the relevant scalar product, in the context of the search for stationary energy levels 
in a static metric. The link with the procedure adopted in Refs. ^Hl ^1 is made 
through the recent interpretation of this procedure by Leclerc [20] • In Sect. El the 
weak-field and non-relativistic limit of the stationary DFW equation in a static metric 
is being studied, i.e., a slow Dirac particle is considered, in a static metric that differs 
little from a fiat metric. The application to ultra-cold neutrons in the Earth's gravity 
is discussed. Our conclusion makes Sect. El 

2 Hermitian Hamiltonian for the Dirac- Fock-Weyl 
equation 

2.1 Starting Hamiltonian for the DFW equation 

The DFW equation may be written as [TH | I2 H 1^: 

(^7^D^ - M) ^ = 0, D^ = d^-T^, M = mc/h, (1) 

where m is the rest mass of the quantum particle, and where the "deformed" Dirac 
matrices 7^ satisfy the anticommutation relation 

7V + = ^g'"' I4, I4 = diag(l, 1, 1, 1), (2) 

with {g'^") = {g^vY^ , 9 = (gfiu) being the (curved) space-time metric, and where the 
matrices are given by 

^, = ^^?A.<,a^s^'^-^^,.,s'^ (3) 

in which the Tx^fj^ 's are the first-kind Christoffel symbols, s^'^ = | (7'*'7'^ — 7'^7^) , and 
the matrices A = (a^) and B = {b^) (with a'j^b^ = 5") transform the natural basis 
= gf^ into the local "tetrad" Uf^ and conversely: up = a^Ca, = b^^Uf^, the tetrad 
being orthonormal: ^ 

9iuf3,u^) = a'^a'^gaf^ = r]f^^, 77 = diag(l, -1, -1, -1). (4) 

^ One may go from the (1,-1,— 1,-1) signature, used here and in Refs. jl5l 1161 FTTl I18^ ■ to the 
(—1,1,1,1) signature, used in Refs. 01l221i by setting: r/' ~ rj, g' ~ —g, and j'^ — The 
bases and the matrices A,B stay unchanged, as well as Eqs. l(5J-(0Ji (of course with primes), while 
Eq. (Q) takes the form assumed in Refs. |22E21- 
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The "deformed" Dirac matrices are related to the "flat" ones 7^^, that satisfy Eq. (0) 
with g = T], by 

7" = 7^- (5) 

Multiplying (P) by 7^^ on the left and using (j2I), one gets the DFW equation in 
Schrodinger form: 

th^ = HV^, (6) 

with 

H-^o'°-'^^Vz), + .fcr, (7) 

(Spatial indices are denoted by Latin letters, space-time indices by Greek ones; and 
a;° = ct.) The Hamiltonian operator ((Tj) makes sense whether Eq. is satisfied or 
not, and may be rewritten as 

H = ». + ^^°-'^^Vfl„. (8) 

The DFW equation (0) is covariant under a general coordinate change: x'^ = (f)^{{x'^)), 
under which the spinor is left invariant |15| I21j. while the 7^^ matrices transform like 
a vector, i.e., 7'^ = j'^dx'^ /dx'^ 21^. (This is because the tetrad is left unchanged in a 
coordinate change.) Hence, j'^Df^tp should also be invariant. It is in fact so: a tedious 
calculation allows one indeed to check that [Eq. hence also D^ip, transforms 
like a covector. It follows then from (jHI that the modified spinor field after applying 
the Hamiltonian operator: liip, remains also invariant under purely spatial coordinate 
changes, x'^ = 0-'((x'^)), = or x' = F(x), t' = t. That is, 

(Hy)(t,x') = (HV^)(t,x), [x' = F(x)], (9) 

with in fact il)'{t,x.') = ilj{t,x.). 

When the metric is diagonal: g = diag(a^), one may define the tetrad from the 
matrix 



A = diag(l/y'|a^|). (10) 
Using this and (jH)), we find that the first sum in Eq. ^ vanishes, and that 

= -i^^g/^- (no sum on /i), S^-" ^ ^ (^^ - TT) ■ (H) 

Let us specialize the Hamiltonian ((Tj) to the particular diagonal metric assumed by 
Obukhov ^: 

ao = [r(x)]^ aj = -[Wi^)]', x^(x^) (j = 1,2,3). (12) 
This is a static metric, i.e., it verifies 

gfiu = gfiui^) and goj = 0, (13) 
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whence g^'^ = Qqq, thus g'^^ = V for the metric ()12p. We find from (fTT|) : 

T/. ■ .n.- 

= 7V, (14) 



^~2iy ~ '^^'='2iy ' ^-[oa'j' ^^^^ 

where cr'(/ = 1,2,3) are the Pauh matrices. (Beginning with (fTH|) . we shall use the 
standard set [1711221 of the matrices 7^.) Putting (fT^ and ^T^i in the Hamiltonian ((7j), 
and using again the diagonal tetrad (fTUI) with (fT^ in the definition (0) of the deformed 
Dirac matrices, we get 



(16) 



Using then the standard definitions of and T} , one checks easily that Obukhov's 
starting Hamiltonian {|I7j, Eq. (12)} is recovered: 

H = m^(3V - ih^ + ^W,^ - ihca^'^d,. (17) 

2.2 Stationary energy levels in a static metric 

For a static metric, the time coordinate t in which we have fll3p is unique (up to a scale 
factor) Thus, if we postulate a stationary wave: 

^(t,x) =0(t)a(x), (18) 

this decomposition (|T8|l must use that "static time" t. {For a time-independent but 
non-static metric, in contrast, the decomposition (fT^ would not be well-defined, since it 
might use any of the different time coordinates compatible with the time-independence 
|25j . say t 01 t' = t + /(x).} Moreover, for any time- independent metric, the Hamil- 
tonian operator ^ does not contain time derivatives, and its coefficients do not de- 
pend on time. Hence, as is well-known, substituting the stationary wave (fTH|) into the 
Schrodinger-type equation ©, leads to 

0(t) = Cexp(-iEt//i), (19) 

while the amplitude function a {which may well have several components, as is the 
case here, provided that the time-dependence (pit) is a scalar [26j} is a solution of the 
eigenvalue problem 

Ha = = huja. (20) 

Equations ()18 p - (j^m) define the stationary energy levels associated with the Hamiltonian 
H — provided, of course, that the "energy levels" themselves, i.e., the eigenvalues E of 
problem ()2()|1 . are real: if that is not the case, the time- dependence ()19|1 contains a 
(real) exponential term, hence the solutions p8|l can in no way be qualified stationary. 
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Hermitian operators make the only well-known class of operators for which we can a 
priori state (thus before calculating the spectrum) that the eigenvalues are real. ^ 

2.3 Scalar product 

The Hermitian conjugate of an operator depends, of course, on the (Hermitian) scalar 
product which is used on the domain of this operator. Thus arises the question of the 
scalar product (a; b) with respect to which H may, or may not, be a Hermitian operator. 
That scalar product should not involve derivatives. Otherwise, when calculating (Ha; b), 
higher-order derivatives would arise. Hence, it has necessarily the following general 
form: 

(a; 6) = [ (a(x).6(x)) dV(x), (21) 

where (y-z) is a Hermitian product defined for arrays y and z of four complex numbers, 
and where V is some volume measure defined on the "space" manifold M, i.e., on the 
preferred spatial section defined by the static metric. The scalar product (y.z) used in 
the integration (j2H) might a priori depend on the integration variable, thus on the point 
X G M, say {y.z)x- However, in the case of the DFW equation, the value ipit,:^.) of the 
spinor is invariant under general coordinate changes [hence the amplitude a(x) = a{x) 
of a stationary ip ()18p is invariant under spatial coordinate changes]. As is usually 
phrased, this means that it belongs to some "internal space" independent of x — say, 
simply to C^. Hence, in the present case, the product {y.z) should not depend on x. 
The only natural product is then the canonical scalar product on C*^, 

(y.z) ^ y^*z^ ^ y^z. (22) 

In contrast, if ip were a 4- vector under coordinate transforms, as is in fact the case for 
the alternative gravitational Dirac equation [14J, the use of (j^^ would be incompatible 
with the coordinate-independence of the product ()2H). The coordinate- independence of 
i]21\i demands also that the volume measure V be invariant under the allowed coordinate 
changes, i.e., under the purely spatial coordinate changes (which are those compatible 
with the static character of the metric). But, for a general static metric, there is only 
one such invariant volume measure, whose expression in a given coordinate system is 

dV(x) = ^/Hx) d^x, h = det{hjk), (23) 

where (hjk) is the component matrix of the metric h on M associated with the static 
space-time metric g (with components g^u), i-e., hjk = —Qjk in static-compatible coordi- 
nates. Thus, we find that there is only one coordinate-independent scalar product that 
appears naturally when investigating the Hermitian character of the Hamiltonian asso- 
ciated with the DFW equation for a static metric — namely, the product (jSI}, specified 
by the definitions ((221) and ((23) • ^ 

^ But, of course, there are plenty of non-Hcrmitian operators with aU real eigenvalues — although, 
within the important class of "normal" operators (i.e., those which commute with their adjoint), an 
operator is Hermitian if and only if its spectrum is real |27j . 

^ Tho ensure that ()21|l is well-defined, one must assume square-integrable functions: a and b G 
L2(M,V). 
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2.4 Link with the hterature 



In Refs. 113 123 Uni CZI, a product with the same form (PT|) . and with (imphcitly) 
the same definition ()22|) for the "local" product, is being used. Each of these works 
considers a particular coordinate system, in which the space-time metric is assumed to 
have a special form — e.g., the form (jl2p in Ref. ^7]. Each in their specific coordinate 
system, these authors use the coordinate volume measure d^x, instead of the invariant 
volume measure (j^Hj) . Thus, they use the scalar product 



ia\b)= / (a(x).6(x)) d^x. (24) 

If one interprets this as the basic, starting, scalar product used by these authors, he 
notes that it is coordinate-dependent, and that the Hamiltonian (|7j) is, in general, not 
Hermitian for this product — e.g., the Hamiltonian (fT7|) is not Hermitian for the product 
(j24j) ■ if the metric has the form (jl2j) . Then, the non- unitary transform used by these 
authors: 

^ = T^/-, H = THT~\ (25) 



with 



[and, in particular. 



> 1/4 

-9 ' 
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9 = det{g^,) (26) 



T = W^l'^ (27) 

for the metric (fT^ ]. is interpreted as a modification of the Hamiltonian with the aim 
of getting a Hermitian Hamiltonian ^H]- This interpretation leads to a problem: the 
modified Hamiltonian H defined by Eq. (^3)) is, in general, not Hermitian any more if 
one transports it to another coordinate system and uses the scalar product defined by 
the same Eq. (j24j) taken in the new coordinates 

However, a more relevant interpretation of the works J3 |2H1 ^1 ^| has been 
suggested recently by Leclerc [20] (Sect. 3 in that reference). He notes that the trans- 
formation (j23)-(|2ni) gets the scalar product 

(a II 6) = \ a(x)t7VKx) d^x (28) 



to the scalar product (j24j) . in an isometric way, i.e. 

(a I 6) = (a II h). (29) 

This is indeed easy to check in the case of a diagonal tetrad (fTUj) (implying a diagonal 
metric), for we get then from Q and 



v^7V = T2l4. (30) 
Due to and (j^ . we have also 

(Ha I h) = (Ha || b). (31) 
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Thus, if the operator H is Hermitian for the product ()28|). then so is also H for the 
product 20 , and conversely; in addition, the two operators have the same spec- 
trum. As noted by Leclerc [201; "the hermiticity properties of the operators can be 
directly read off' with the flat product (j^ : for example, see Eq. (14) in Ref. [T7] . 
The scalar product (j28j) arises from the conservation law associated with the standard 
DFW current, and it may be rewritten in a coordinate- independent way (for a 
given hypersurface, defined in the starting coordinate system by = Const.). Here, 
we show that the product coincides, in the case of a diagonal tetrad (fTUj) . with 
the coordinate-independent product (j2T| - (j23|) that we introduced. We have in the most 
general case j2S]: 

T^ = Zl = h. (32) 
fi'oo 

For a diagonal tetrad, we have Eq. (jHUj) . whence indeed from (|^: 

(a II 6) = (a; 6). (33) 

We conclude that, in the case of a static metric which can be reduced to the diagonal 
form (which is the case considered here and in Refs. IT7]). the invariant scalar 
product ()2H) - ()23|1 can be rewritten in the form ()28j) . Then, the non-unitary transform 
fj25p is indeed "a purely mathematical operation which does not affect the physics," 
and which transforms that product to the fiat product ()24|). in terms of which the 
hermiticity of the Hamiltonian operator is easily checked [20] • In particular, in the case 
of the metric p2|l . this transformation brings the Hamiltonian (fTTj) to an explicitly- 
Hermitian form ^^1, hence we know that the starting Hamiltonian (fTTj) is Hermitian 
for the invariant product (PT|) - (f^ . Nevertheless, it is interesting to show that the 
adjoint of the Hamiltonian (fT7j) for that product can be determined directly. ^ 

2.5 Adjoint of the Hamiltonian (direct method) 

Thus, let us compute the adjoint H^ of the operator H, for the scalar product ()2H1 - 
fj23p . To do that, we need to know the adjoint operator of the partial derivative dj 
(or, equivalently, that of the momentum operator pj = —ihdj). Of course, the operator 
dj is coordinate-dependent, and so will be its adjoint. But, since the Hamiltonian ((Tj) 
verifies Eq. its application to a merely space-dependent spinor a(x) verifies 

H'a'(x') = Ha(x). (34) 

The scalar product (|2H ) - (j23j) being invariant too, the resulting adjoint operator H^ is 
also invariant. I.e., if one makes an admissible (thus purely spatial) coordinate change, 

^ When investigating the hermiticity of the Dirac(-Fock-Weyl) Hamiltonian, or, more exactly, when 
checking whether (Hip \\ ip) = {ip \\ Kip), Leclerc [201 assumes that ip verifies the wave equation. 
However, this identity [and in fact the more general identity (Rip \\ ip) — (ip \\ Hi^s)] has to be verified 
by general functions ip (and Lp) in the domain of H, not merely by ones satisfying the wave equation. 
Moreover, in the context of the search for stationary solutions H18|l . the relevant hermiticity applies to 
spatial wave functions a. The associated time-dependent wave function (|18|l obeys the wave equation 
iff a is an eigenfunction of H, hence checking the hermiticity for such wave functions is not enough. 
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x' = F(x), t' = t, then the adjoint operators will correspond by Eq. (jH^ (with in the 
place of H). We determine first the adjoint dj of the operator dj, considered as acting 
on scalar functions, thus with respect to the invariant scalar product that corresponds 
to (PT|) for scalar functions, 

(a, b)= I a*h dV. (35) 



M 



If the scalar function / cancels outside a bounded domain, or more generally if [in the 
coordinates {x^)], f{^)\Jh{^ = o(l/r^) as r = |x| — > cx), ^ then one has 

L Try '^=Lii [f(-)VH^] = 0- (36) 

(This is got by applying the divergence theorem to the vector field /(x)a//i(x) e^- with 
the ball B(0, r) in the Euclidean space R^.) By using this result with / = a*b, we obtain 



M 



(a*,6 + a*b ,) dV + /" a*6^ dV = 0. 



(37) 



As we know, the adjoint of an operator A [here one acting on scalar functions, thus 
with the scalar product ()35|) ] is defined to be the operator A-^- such that 

(Va) (V6) {a,Ab) = {A^a,b). (38) 

Hence, it follows from Eq. ()37|) that ^ 

S] = -3, - |. (39) 

with, for the metric ()12|l . 

hi Wi , , 

2f = V- («) 
It is easy to check that, when an operator A is extended from scalar functions to 
"vector" ones by A.(a^) = (Aa^), its adjoint A''" for the product (f7T|l - (j^ is obtained 
by extending in the same way the adjoint A^ for the product Also as expected, 

the adjoint, for the product (PT|) - (j^^ . of the operator defined by a complex matrix M, 
is the operator defined by the adjoint matrix, M'^ = (M*)""". Thus, in Eq. ()17|) . the first 
operator is Hermitian, the second one is anti-Hermitian, and the adjoint of the third 
one is [Eqs. dHHIl-dini)]: ^ 



a-' 



<Fd,+F,)-?>z^F 



(41) 



^ Thus, we assume that the space M is diffeomorphic to R^, and that the chart x ^ [x^) covers 

M. 

^ Without going into detailed considerations on the domains of the operators, it is clear that the 
functions to which H37(l applies will be dense in the relevant domain. 

^ The usual convention is followed for product operators: thus, the product operator djF transforms 
a to dj{Fa) = {Fdj + Fj)a, with Fj the mere multiplication by the function dF/dx^ . 
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We find then easily that 



Ht = mc^pV - ihca^ [Fdj + Fj + 2^F - ^ j (42) 



"->--«-'(^ + ^ + ^9.)=H. (43) 



Thus, we checked directly that the Hamiltonian operator ()17|) is Hermitian for the 
invariant scalar product ()2H) - (j^^ . Therefore, to compute the stationary energy levels 
of the Dirac particle, we may consider the eigenvalue problem (j2Up with this starting 
Hamiltonian (fTTj). Using the standard set [T71I2S1 of the matrices and /3, and writing 
the 4-component amplitude function a(x) as a couple of 2-component functions, 

a(x) = (^(x),x(x)), (44) 

we may write explicitly ()20|) (after division by he) as 



Thus, X is expressed as a function of by Eq. and is eliminated by reporting 

this in ()45|1 . In this sense, the stationary DFW equation also has only two independent 
complex degrees of freedom, as is well-known for the flat-space-time Dirac equation. 



3 Post-Newtonian approximation (PNA) for station- 
ary energy levels 

3.1 Framework: asymptotic PNA 

A PNA was already considered in the quantum domain by Kiefer & Singh 29 , who 
expanded in powers of the phase function in the "minimally-coupled" gravitational 
Klein-Gordon equation, in a non- stationary situation. This approach was further de- 
veloped by Lammerzahl jSU], who applied it to explore the modification of the coupling 
of matter to the electromagnetic field, which is induced by the presence of a (weak) 
gravitational field. Here, we are using the so-called "asymptotic" PNA [^ESESl, i-e., 
we consider a (conceptual) family (Sa) of systems. The aim of this approach is to give a 
mathematically clearer meaning to the postulated expansions, interpreting them indeed 
as asymptotic expansions. In the present case, each system Sa is constituted by a mas- 
sive body producing the static metric, and by a Dirac particle which is in a stationary 
state in this metric. Since we are studying a "test quantum-particle", assumed to not 
affect the metric, we need to know only the metric produced by the body: let g^^^ be the 
metric in system Sa. We do not need to specify which theory of gravitation is assumed. 
We suppose that there is a coordinate system t,x, such that, for any A, the metric g*^^^ 
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has the form p2|) with the following (post-Newtonian) asymptotic expansion as A — * 0: 

yW(x)^ Ujf^' = l-A^ + 0(A^), (47) 
L J 

lyW(x) = [-n^f 1 = 1 + A + 0(A2) (no sum on j). (48) 



c2 



This will indeed be the case if a post-Newtonian family of static bodies is envisaged 
either in the framework of GR in the harmonic gauge [S3innilSlE2]! or in the framework 
of a recent scalar theory . In both cases, U is formally identical to the Newtonian 
potential: it is defined in the coordinates t, x by 

f/(x) = G J p(y)dV|x-y|, (49) 

where p arises as the first coefficient in the expansion of the energy density T°°*^^) in 
the coordinates t, x, with T^"**) the energy-momentum tensor in system S^: 

T^oix) = Ap [1 + 0(A)] . (50) 

Note that, therefore, the first approximation of the energy density in system Sa (the 
"Newtonian density") is p\^^ = Ap. Hence, the real equivalent of the Newtonian poten- 
tial is not U (which does not depend on A), but is instead, in system Sa: 



4')(x) = At/(x) = p^")(y) dV/ |x - y| 



(51) 



The system of interest, S, e.g. the Earth (assumed isolated), is supposed to 
correspond to a small finite value Aq -C 1 of the parameter: S = Saq, so that it 
makes sense to use the asymptotic expansions for that system. This means that the 
gravitational field in system S is indeed a weak field in the physical sense. ^ Moreover, 
the classical velocity of the test particle is assumed to be of the same order in the small 
parameter as that of a test particle orbiting in the weak gravitational field considered, 
namely |H11 EH EH El 

u(^)(t) =ord(v^) = uo(t)v^ + 0(A3/'). (52) 

This means physically (in system S) that the initial velocity of the particle (e.g., the 
initial velocity of the neutron fiux) is at most of the order of orbital velocities (as is 

^ This is not automatical if one considers an a priori given, abstract family (Sa), as is done in 
Ref. because asymptotic expansions like (|47|I - H48|I remain true if one changes A to A' = a\, 

with a a constant: in that case, the fact that A is small for some particular system S — Saq has 
no objective meaning. But in physical practice, one defines Aq in the very system of interest S, e.g. 
Aq = SuPxgR3(l — V{x.)) for the metric ((T^ . One must check that it is indeed negligible w.r.t. the 
unity: this means really that the gravitational field in system S is a weak field in the physical sense. 
Then, the family (Sa) is deduced from the data of S [62i (This is done by deducing a family of 
initial data from the initial data corresponding to S.) 
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certainly the case in the experiments [3131^1); so that it makes sense to use the PNA — 
then, u (and also u) will remain of PN magnitude as the time goes. It follows that the 
classical energy of the particle in the gravitational field [23] admits the expansion 

41.a. - {giS + {P'''?c'] = mc' + 0(A). (53) 

3.2 Application to the stationary energy levels 

The stationary energy levels of the particle, considered as a Dirac (quantum) particle, 
must have the same order of magnitude as the particle's classical energy. Hence, any 
solution E^^^ of the eigenvalue problem (jSDllA (i-e., for system Sa), should have an 
expansion 

= mc^ + 0(A) (54) 

(assuming that one may thus follow as a function of A each of the different eigenvalues). 
This asymptotic expansion of an energy level implies that it is positive (as soon as A 
is small enough). More precisely than (jH^ . we shall assume a first-order expansion for 
E^^\ or rather equivalently for u^^^ = E^'^^ /h: 



c 



M + cuiA + 0(A2), (55) 



as well as for the amplitude function a^^^^ = {^p^^\ x^^'^)- 

(^W = ^0 + ^iX + OiX'), x^'^ =Xo + XiA + 0{\'). (56) 

Inserting and into the explicit eigenvalue equations and (jlHI), using the 
expansions (fTTj) and (PHj) of V and W, and identifying powers of A, yields, at the order 
zero, 

-tcr^Xo,j = 0, (57) 

-zaVo,, = 2Mxo, (58) 
and at the order one (with U = U/(?): 

-%a' (xi,o - 2UXo,j) - \(y'U,jXo = (^1 + M[/) (^o, (59) 



-la^ ((^i,,- - 2f/(^o,i) - 'i^^'U^m = 2Mxi + (^1 - Mt7) xo- (60) 
Using the well-known property of the Pauli matrices: 

= ( f ^: ^ f (61) 

[ I2 if J = /c, 
we note first that (j57j) and (j58|) imply 

A(^o = 0. (62) 
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Then we eliminate xo and xi from ()59|) by using 
help of dSZI), dnH) and 



and (jUUj) . and we get with the 



2M 



3z ~ 1 ~ 



oJi + MU]ipo. 



(63) 



In Eqs. (jH^ and ()63|1 . we have four independent scalar equations for the four scalar 
unknowns contained in the zero-order and first-order coefficients ipo and ipi of the expan- 
sion of the 2-component field y?. This makes the IPN eigenvalue problem determinate. 
However, we may use Eq. (jU^ so as to rewrite (jU^ as a single approximate equation 
for the IPN field (which is the order-one approximation of v^''^-*) 



V^(i) =Vo + Av^i- 
Owing to and (jBH), we have indeed in Eq. 

In the same way, it follows from ()55p and ()64j) that 



(A) 



M)y,|^; + 0(A). 



(64) 



(65) 



(66) 



Thus, multiplying (jH!^ by \hc^ we get [omitting the superscript ^^"^ and the subscript 
(1), i.e., defining ({) 



<^|JJ, and also = U^^'^ = XU, Eq. 



- 7^^'^ - mUf^if + - — - 
2m Zmc^ 



3 11' 
(E-mc2)v9 + 0(A2). 



This equation shows explicitly the additional terms with respect to the stationary 
relativistic Schrodinger equation in the Newtonian gravity field. 



-—Ay. 



(67) 



non- 



(6^ 



(The subscript "nr" stands for "non-relativistic") It is interesting to note that, in Eq. 
(jnZj), all terms are of the same order in the small parameter A. [This is a priori obvious 
since ()67|) is a mere rewriting of the exact equation ()63p . the latter involving only 
expansion coefficients like y^o, v^i and f/, which are, of course, of order zero in the small 
parameter.] Namely, all terms are order A. ^ This result means that, from the point of 
view of the asymptotic FN scheme, the corrections to the non-relativistic Schrodinger 
equation (fH^ — e.g., the corrections to the energy levels — do not need to be small with 
respect to the corresponding quantities as they are found using Eq. (|UHj). However, an 
asymptotic estimate is not a numerical one: we must now numerically investigate the 
differences in the energy spectrum. 

^ On the r.h.s. of H67|l . we consider {E — ■mc^)Lp as one term: it is of order one in A, but it is the 
difference between two terms of order zero in A. 



13 



3.3 Estimates for ultra-cold neutrons in the Earth's gravity 

The mass of the neutron is m ~ i x 10"^^ kg (mc^ ~ 939.57MeV ~ 1.50 x 10"^° J). 
The Earth's Newtonian potential is, on Earth, ~ ^ ~ ^'^^''l^jxw'^''^^'^ - 
6.3 X 10''m^/s^, whence mll^ ~ 10""^^ J. We have |Vt/N| = 9 — 9.81 m.s~^. Assuming 
perturbatively for ip the functions found from the non-relativistic Schrodinger equation 
for the first energy levels of ultra-cold neutrons in the Earth's gravity ^6^, we have 
iVv^l ~ 10^ m~^. Now the relativistic corrections involve the minute coefficient 

(1.054 x 10-34)2x3x1026 .q, 2 

C = ^ ~ ^ '-Ta 3.7 X 10-^^ kg.ml 69 

^ 2mc2 9 X 1016 ^ ^ ^ 

Therefore, the two first correction terms in the square bracket in Eq. (p7j) . which are 
of the order of C,g | Vv?| — 3 x lO^^^J, are utterly negligible with respect to mU-^ip ~ 
10"^^ J. It is even more so for the third, last correction term in (pTTjl . —^ipAU^ with 
-Af/N = 47rGpN ^ 8 X lO-^s'^ (with ^ lO^kg/m^). 

However, one should account for the fact that the energy levels are defined only 
up to a constant: in Eq. (j67j) . one may change simultaneously to = + C and 
E to E' = E — mC. [The same is true in Eq. (jHEI)-] It follows that, in the experiments 
with neutrons passing through a horizontal slit [Hj, which are very localized as compared 
with the Earth's radius, we may approximate the Newtonian potential as 

f/N ^ -gz, (70) 

with z the vertical coordinate, counted upwards from the lower side of the slit, and 
taking values up to a few 10"^ m. This gives now —mU^ ~ 10"'^'^ J. That has the 
same magnitude as the non-relativistic energies E^j. E ~ mc^ [S]. It still exceeds the 
relativistic corrections by some 23 orders of magnitude. 



4 Conclusion 

In this work, we derived the post-Newtonian equation for the stationary energy lev- 
els of a slow Dirac(-Fock-Weyl) (DFW) particle in a weak static gravitational field, 
Eq. (jnZj). To our knowledge, this equation was not derived before, although weak-field 
expressions for the DFW Hamiltonian have been proposed fSJ El El QB] • In a first 
step, we showed that there is just one natural coordinate-independent scalar product 
relevant to stationary wave functions in a static metric, given by Eqs. (PT|) - (j^^ . In 
the case of a diagonal tetrad, this scalar product coincides with the scalar product (PH|) 
considered by Leclerc [201 , and which arises from the conservation law of the DFW 
current. We checked directly that, at least in the case where the metric can be set in 
the "space-isotropic" form ()12|) assumed by Obukhov T7 , the Hamiltonian ((7j) of the 
DFW equation turns out to be Hermitian for that scalar product. Instead, de Oliveira 
& Tiomno (T^j, Varjii & Ryder [12] and Obukhov [T7j used the non-unitary transfor- 
mation (j2Sl)- That transformation gets the scalar product ^2E\i to the "fiat product" 
(|24j) . with which the hermiticity of the modified Hamiltonian, equivalent to that of the 
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starting Hamiltonian for the product ()28|). is easily checkable. 

In a second step, we used the "asymptotic" scheme of post-Newtonian approxima- 
tion (PNA) jSIlE2lE3] to determine the approximate equation that governs the station- 
ary energy levels, in order to be able to compare that equation with the non-relativistic 
Schrodinger equation in the gravity potential. To use that scheme, we assumed that the 
gravitational field is weak and that the classical velocity of the Dirac particle is small, 
both in an asymptotic sense, i.e., for A —>■ — considering a family (Sa) of systems. 
It makes sense to apply the results to the system of interest (e.g., a flux of ultra-cold 
neutrons in a laboratory on the Earth), insofar as the relevant value Aq of the param- 
eter for that system is small (as is indeed the case in the example mentioned). In our 
final equation ()67p. the new terms, as compared with the non-relativistic Schrodinger 
equation (jUHj) . are of the same order in the small parameter as the terms which are 
already there in Eq. ()68|) — namely, they are order one. This is a surprising result if 
one compares it with the situation for classical particles, in which situation the first 
relativistic corrections are one order higher in the small parameter than the Newtonian 
terms. However, this result does not mean that the corrections have the same numerical 
magnitude as the "main" terms. And indeed, we found that, for ultra-cold neutrons in 
the Earth's gravity field, the corrections of the DFW equation to the energy levels are 
hopelessly negligible with respect to the energy levels as found with the non-relativistic 
Schrodinger equation. 

We note that, in their works on the nonrelativistic limit of the DFW equation, 
de Oliveira & Tiomno jT3], Varjii & Ryder [12], Obukhov [T7], as well as Silenko & 
Teryaev [TH] , used Foldy-Wouthuysen EE] transformations, or at least (according to 
the discussion in Ref. [18j), transformations of the Foldy-Wouthuysen type. In contrast 
with these authors, we consider explicitly the stationary energy levels, and this in the 
post- Newtonian approximation. These two points, taken together, allow us to state Eq. 

and to use it, thus automatically selecting the positive-energy solutions. For this 
reason, we did not need to have recourse to a Foldy-Wouthuysen transformation. 

Acknowledgement. I am grateful to F. Selleri, F. Romano, to the late C. De 
Marzo and to his nice team, for their warm hospitality in Bari. Also, thanks are due to 
Franco Selleri and to Konstantin Protasov, who urged me to evaluate the magnitude of 
the corrections to the energy levels. Finally, the referee pointed out the very relevant 
work 120] and asked questions about the change of units used in the previous version 
|19j . which led me to realize that this change is unnecessary in the present case — thus 
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